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We revisit the issue of color neutrality in effective model descriptions of dense quark matter 
based on global color symmetry. While the equilibrium thermodynamics of such models is now well 
understood, we examine the collective modes, focusing on the fluctuations of the order parameter. 
We point out that the constraint of color neutrality must be carefully generalized in order to obtain 
physically consistent and well-defined results. Particularly important is that the collective modes 
associated with order parameter fluctuations couple to charge density fluctuations in the neutral 
medium. We start by proving explicitly that, in contrast to claims made previously in literature, 
Nambu-Goldstone bosons of spontaneously broken global color symmetry remain exactly massless 
even after imposing the color neutrality constraint. As the next step, we make the argument general 
by using effective field theory. We then employ the high-density approximation to calculate the 
couplings in the effective Lagrangian and thus the Nambu-Goldstone boson dispersion relations. 
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I. INTRODUCTION 

The physics of cold dense quark matter is governed 
by the theory of strong interactions, the Quantum Chro- 
modynamics (QCD). Unfortunately, it becomes strongly 
coupled in the phenomenologically interesting range of 
densities. At present, there is no analytic method in 
the market that would be able to perform reliable first- 
principle calculations in this regime, although first at- 
tempts have already been made [l|. Likewise, lattice sim- 
ulations at high density are out of reach of the standard 
Monte-Carlo techniques due to the sign problem. 

As a consequence, one usually has to resort to sim- 
plified calculations using models which more or less imi- 
tate the full QCD. Since the technical difficulty of QCD 
stems from the strong gauge interaction, these models in- 
variably replace the color gauge symmetry with a global 
one, the prototype being the Nambu- Jona-Lasinio (NJL) 
model (see Ref. [2| for extensive reviews). One then 
has to deal with various model artifacts that are par- 
ticularly severe at high baryon density where cold quark 
matter is expected to behave as a color superconductor 
(see Ref. [3| for recent reviews). A condensate of quark 
Cooper pairs breaks color symmetry and (in most color- 
superconducting phases) induces nonzero color charge 
density. This cannot be physical since, due to the long- 
range nature of gauge interactions, it would give rise to 
nonextensive energy density and an ill-defined thermody- 
namic limit [4]. In QCD, it is compensated by an induced 
gluon condensate so that the system as a whole is color 
neutral @- On the other hand, in models with global 
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color symmetry color neutrality has to be imposed as a 
thermodynamic constraint. 

Arranging for global color neutrality by introducing 
one or more chemical potentials associated with the color 
charge(s) costs energy as compared to the unconstrained 
equilibrium state. The neutrality constraint then has to 
be defined carefully in order to avoid spurious instabili- 
ties. Initially, only neutrality with respect to the three 
colors of fundamental quarks was required in literature 
on color superconductivity, leading to the existence of 
seemingly neutral and energetically preferred states [6(. 
As was pointed out in Ref. [7| , such a restricted neutral- 
ity requirement is sufficient only for special orientations 
of the diquark condensate in the color space. In general, 
neutrality with respect to all eight generators of color 
SU(3) has to be imposed, and the full set of eight color 
chemical potentials then have to be introduced Q. 

Once the chemical potentials are fixed to make the 
equilibrium ("ground state") neutral, color symmetry is 
apparently broken explicitly. This observation led to the 
conclusion that the Nambu-Goldstone (NG) bosons of 
the spontaneously broken global color symmetry have 
small, yet nonzero masses, proportional to the color 
chemical potential(s) Q . The goal of the present paper is 
to show that this conclusion is premature. The physical 
picture behind our claim is as follows. The introduction 
of chemical potentials is enforced by the diquark con- 
densate, which itself is generated dynamically as a non- 
perturbative solution to the equations of motion. The 
color chemical potentials therefore are not mere external 
fields, but stem from the dynamics of the system. In full 
QCD, this is indeed the case [5]. After all, none of the 
color symmetry is broken explicitly once the dynamics in 
the gauge sector is properly taken into account. In the 
NJL model, the color chemical potentials mimic the role 
of the gluon condensate in the full QCD. 



One should note that the chemical potentials com- 
monly used in literature are demanded to make the 
ground state neutral. However, the NG bosons constitute 
its excitations. One therefore cannot use the very same 
values of the chemical potentials when dealing with such 
non-equilibrium field configurations, and the thermody- 
namic constraint of color neutrality needs to be extended 
properly. We show that this can be technically achieved 
by treating the chemical potentials as dynamical vari- 
ables, or in other words, secondary order parameters in 
the gauge sector induced by the primary order parame- 
ter, that is, the diquark condensate. Once this is done, 
the NG bosons remain exactly massless, as predicted by 
the Goldstone theorem. 

The plan of the paper is as follows. In Section [TT| we 
explain the central idea of the paper using a very sim- 
ple toy model. While this is rather trivial, it is intended 
to demonstrate the conceptual simplicity of our strat- 
egy, which might otherwise be concealed by unimportant 
technical details of the NJL calculation to follow. Sec- 
tion IIIII provides a correction of the calculation of the 
NG boson masses of Ref. [9( , showing that they are ex- 
actly zero once color neutrality is properly imposed. In 
Section llVI we extend the calculation and determine the 
full dispersion relations of the NG bosons. We employ the 
high-density approximation [k| (see Ref. [11| for a review 
and further references), which both simplifies the calcu- 
lation and makes the results model independent, for it is 
known to capture the leading order of the high-density 
asymptotic behavior in full QCD. Finally, in Section W\ 
we summarize and conclude. 



II. NEUTRALITY IN A SCALAR TOY MODEL 

Consider the scalar theory with a global SU(2) x U(l) 
symmetry, defined by the (Minkowski space) Lagrangian 



C = V^Vcj) - M- 



A 



(1) 



where ^ is a complex doublet field. This model was 
investigated many times before [12| and it was shown 
that when the symmetry is spontaneously broken in pres- 
ence of a chemical potential associated with the U(l) 
subgroup, the three broken generators give rise to one 
type-I and one type-II NG boson whose dispersion re- 
lations at low momentum are linear and quadratic, re- 
spectively |13| . This is in accordance with the general 
Nielsen-Chadha counting rule for the number of NG 
bosons as well as the fact that in the ground state, isospin 
acquires nonzero density [14J . 

Here we wish to investigate the effect of enforcing neu- 
trality with respect to the isospin SU(2) group. This 
may be regarded as a simple toy model for understanding 
charge neutrality in non-Abelian gauge theories where 
the global symmetry (after gauge fixing) is spontaneously 
broken such as in color superconductors. We will explain 
that using the values of chemical potentials obtained from 



the neutrality constraint on the ground state to deter- 
mine the excitation spectrum leads to a spurious insta- 
bility. We will then use two different approaches to deal 
with it. First, we shall demand, and justify, that all 
uniform field configurations be SU(2) neutral, whether 
they correspond to the thermodynamic equilibrium or 
not. A more model-independent argument will be given 
afterwards, showing that the SU(2) chemical potentials 
can be treated as induced, secondary order parameters. 
These two approaches will be later used respectively in 
Sections IIIII and IIVI Here we just remark that while the 
former has the advantage of being conceptually more 
straightforward, the latter is more powerful and allows 
us in particular to determine the full dispersion relations 
of the NG bosons. 



A. Spurious instability and its cure 

We start by introducing independent chemical poten- 
tials for all generators of the symmetry group, /x for U(l) 
and jJ for SU(2). The covariant derivative in Eq. ([1]) then 
reads 



V^cp = [<9 M - iS^ (n + f ■ p)]<t>, 



(2) 



where f is the vector of Pauli matrices. The presence 
of chemical potentials for all charges allows us to easily 
find the associated charge density operators (here for the 
isospin charges), 

™=^= l {^Td 4> - do^T(t>) + 2^(^ + p) ( f>, (3) 

Using this expression, the Lagrangian ([T]) becomes 

C =d^d> l <l> + ifi(tfd (t> - d (f> f (f>) + fi ■ n 
-tf(M 2 -fj?+jf)(f>-\( 



(4) 



All matrix structure is now hidden in the density n. 

Let us assume that the scalar field develops nonzero 
vacuum expectation value (the conditions for this to hap- 
pen will be discussed below) and choose the ground state 
as usual as 0o = (0, v) T . The complex doublet field will 
thus be parameterized as 4>\ — ip, <j>2 — v + H + id, 
where ip is a complex field, whereas H and 9 are real. 
For constant field configurations, the charge densities © 
are then expressed as 



—ri\ — /ii(w 2 + 2vH) + 2/iiiRey + second-order terms, 
—rT-2 = fi>2(v + 2vH) — 2fjLv\nnp + second-order terms, 



2«3 = (Ms - V)(v 2 + 2vH + H 2 + 9 2 ) + (/i 3 + /i W 



(5) 



Observe that when ft = 0, the isospin density, n^ = 
— 2/iu 2 , appears as soon as the field condenses. For 



M 2 > 0, this happens once the chemical potential [i ex- 
ceeds the mass M. This is the Bose-Einstein condensa- 
tion. 

If we now want to make the ground state SU(2)- 
neutral, the first terms in the expressions for charge den- 
sities lead to the conditions 



Mi = A*2 = 0, 



/13 = [i (neutral ground state). (6) 



The static bilinear (mass) part of the Lagrangian §4§ 
reads, up to an overall minus sign, 

U bilin (0) =2v(M 2 + 2Xv 2 )H + (M 2 + 6Xv 2 )H 2 (7) 

+ (M 2 + 2Xv 2 )9 2 + (M 2 + 2Xv 2 - V^V 

From the linear term we see that when charge neutrality 
is imposed, spontaneous symmetry breaking can only oc- 
cur for M 2 < 0, that is, when it already appears in the 
vacuum. This is natural: normally, Bose-Einstein con- 
densation (at zero temperature) sets when charge density 
starts to be nonzero, but here we keep the isospin density 
equal to zero, which prevents the field from condensing. 

Substituting the vacuum expectation value, v 2 — 
—M 2 /2X, we find that the Higgs mode has the mass term 
AXv 2 H 2 , while the phase fluctuation of the condensate, 
9, becomes massless. However, the same is not true for 
ip. This observation led in the context of color supercon- 
ductivity to the conclusion that the chemical potential 
needed to render the so-called 2SC phase color-neutral 
breaks some of the generators explicitly and the associ- 
ated NG bosons thus naturally acquire nonzero mass |9j . 
In our case this corresponds to the /13 chemical poten- 
tial "explicitly breaking" the 7^2 generators, i.e., giving 
mass to the tp excitation. However, the situation is even 
worse: the mass squared of ip is negative! A similar prob- 
lem appears in the 2SC phase, as will be shown in the 
following section. In fact, the presence of terms in the 
Lagrangian Q with a single time derivative makes the 
discussion slightly more complicated than just concluding 
that ip has a negative mass squared. One finds that the 
(anti)particle mode annihilated by ip (ip^) has the disper- 
sion Ej- = |fc|=p2/i. Consequently, the particle mode with 
momentum smaller than 2/z is "tachyonic" . To conclude, 
fixing the chemical potential to make the ground state 
neutral as we just did obviously leads to an unphysical 
result. 

The roots of the problem lie in the fact that we fixed 
the chemical potentials in the ground state once for all. 
The fluctuations of the order parameter then drive the 
system off the neutrality, thus naturally lowering the en- 
ergy. This feature was already observed in Ref. [6| , where 
it was erroneously interpreted as an instability of the 
2SC ground state. Buballa and Shovkovy [7]] pointed 
out that the instability disappears when the chemical po- 
tentials necessary to make the ground state neutral are 
transformed simultaneously with the ground state itself. 
It should be stressed that our approach is technically 
very close to that of Ref. [7J : after all, the NG collective 
modes in the infinite- wavelength limit correspond merely 



to a change of the ground state. However, we make one 
step further by demanding charge neutrality also for cer- 
tain non-equilibrium field configurations, corresponding 
to collective modes in the infinite- wavelength limit. In 
other words, we will require that for every uniform field 
configuration the chemical potentials acquire such values 
that the system has zero SU(2) charge. This is reason- 
able since otherwise there would be a uniform charge dis- 
tribution, leading to an ill-defined thermodynamic limit. 
Such considerations are sufficient to determine the mass 
spectrum of the theory, and we will now demonstrate 
that the NG bosons are indeed exactly massless as they 
should. A formalism how to deal with nonuniform field 
fluctuations, that allows us to determine the NG boson 
dispersion relations, will be developed below. 

Let us therefore assume constant fields H, 6, <p and de- 
mand that the SU(2) charge densities ([3]) are zero for all 
values of the fields. This gives the chemical potentials 
the following values, 



2jx 2// 

Mi = He ip, fi-2 = -\ Im 93, 

v v 

Us = u off (neutral excitations), 

v 



(8) 



to lowest nontrivial order in the fields. Substituting this 
back into Eq. (j4|) , the term pi ■ n vanishes by definition of 
charge neutrality, while the rest reduces to 



Vbiiin(0) =2v(M 2 + 2Xv 2 )H + (M 2 + 6Xv 2 )H 2 
+ (M 2 + 2Xv 2 )6 2 + (M 2 + 2Xv 2 )iptip. 



(9) 



Upon minimizing the potential to make the linear term 
vanish, the field ip> becomes exactly massless as it should. 



B. Effective field theory approach 

The masslessness of both 9 and ip can be achieved nat- 
urally once the chemical potentials are treated as dynam- 
ical fields and the local symmetry transformation which 
generates the NG excitations is accompanied by a corre- 
sponding transformation of the chemical potentials. The 
validity of the Goldstone theorem is then an immediate 
consequence of the exact symmetry of the action. Let 
us explain on our simple example the line of reasoning. 
In the underlying gauge theory the charge neutrality is 
maintained by a gluon condensate which compensates for 
the charge of the primary order parameter, here the vac- 
uum expectation value <po- Under a symmetry transfor- 
mation, the primary order parameter and the gluon con- 
densate transform simultaneously. In order to correctly 
capture the symmetry properties of the gauge theory in a 
model with a global symmetry, we must allow the chem- 
ical potentials to transform in the same way the gluon 
condensate would. 

Technically this means that we deal with a theory 
with a global symmetry and two order parameters: the 



primary <po, and the secondary one, that is, the chem- 
ical potentials. A low-energy effective Lagrangian for 
the NG bosons is then constructed as usual by per- 
forming a spacetime-dependent symmetry transforma- 
tion on the order parameter(s). In our specific exam- 
ple, let us parameterize the scalar field as 4> = U(n)(j)o, 
where U(tt) = exp(-r • 7?). Analogously, the tempo- 
ral background gauge field A that contains the chem- 
ical potentials and enters via the covariant derivative, 
T>fj_(/) = (9 M — i8^qA)4>, is parameterized as 



A = U(ir)A U{ir) 



(10) 



where Ao = \x + r • fi with the ground state values of the 
chemical potentials determined by Eq. (J6)). Note that 
this parameterization abandons the amplitude mode H , 
that is, it defines the nonlinear sigma model. 

Inserting the parameterization into the Lagrangian ([1} , 
one finds 



- i<t>l(d U 1 U)A (j>o + <plA Q A Q (t) 
o-A(^ ) 2 . 



(11) 



AP 



It is easy to see that the second, third, and fourth term 
drops because Ao4>o = 0, while the last two terms do 
not include the NG fields 7?. The Lagrangian therefore 
reduces to C = ^(d^U* 1 )(d tl U)(f>o up to a constant. Inter- 
estingly, it does not depend on the chemical potentials at 
all. In particular we can see that all NG bosons if are ex- 
actly massless and have the usual Lorentz-invariant dis- 
persion relation. This is due to the fact that despite the 
chemical potentials, the densities of all conserved charges, 
including the U(l) charge with respect to which neutral- 
ity is not required, are zero in the ground state |15| . 



III. NG BOSON MASSES IN NJL MODEL 

The aim of this section is to show that when color neu- 
trality is implemented using the strategy outlined above, 
the NG bosons become exactly massless in the NJL model 
description of color-superconducting quark matter. Our 
plan is as follows. We first introduce the model and es- 
sentially repeat the calculation of Ref. [9] for fixed color 
chemical potentials. Next, we adapt the approach intro- 
duced in Section lll Al for the present purposes. Finally, in 
Section IIIIC I we present some details of the computation 
of the corrected mass spectrum of the NG bosons. 

For the sake of simplicity, we will use a version of the 
model considered in Ref. [9] which does not take into ac- 
count dynamically generated constituent quark masses. 
This cannot change the conclusions qualitatively. More- 
over, in color superconductors the constituent masses of 
u and d quarks are usually small. Since the calculation is 
rather technical, we will omit details, yet providing the 
key steps and definitions necessary for the reader who 
might desire to reproduce our results. 



The model is defined by the Lagrangian in Minkowski 
space, 

£ = q(i0 + lo VaT a -m)q+ -J2(fP a q)(qP a q C ), (12) 



where q c = Cq T is the charge conjugated Dirac spinor 
and the matrices P a specify the structure of Cooper pairs. 
In the 2SC phase, they read (P a Yb c = l^taba where 
i,j are flavor and a,b,c fundamental color indices, re- 
spectively. Also, P a = 7o-Pa7o- We introduced nine 
chemical potentials fi a , a = 0, . . . , 8, associated with the 
symmetry generators T a in the color space, normalized 
by tr (T a T{)) — ^S a p. They include the quark number 

chemical potential \x, related to \lq by \x§ = /i-^/6 since 
To = 1/V6. 

In the mean-field approximation, the thermodynamic 
potential Q of the model is given by 

(3n = pV^±-±Tr\ogS~\ (13) 

where V is the space volume and A Q is the collective 
field that stands for the composite operator tfq P a q. Its 
ground state expectation value is determined by the min- 
imization of the thermodynamic potential. Furthermore, 
the quark propagator in the Nambu space, \]/ = (q, q c ) T , 
reads 



s-Hk) 



lo^ a T a 
A*P 



&aPa 

- 7oMQ r J - m 



(14) 



The symbol "Tr" in Eq. (jT3|) denotes a trace in the op- 
erator sense. 

Let us introduce some further notation. In the follow- 
ing, we will not need to work with the general orientation 
of the condensate A Q in the color space as well as with 
all the chemical potentials [i a . In fact, all integrals will 
be evaluated with only the A3 = A and /Uo,8 compo- 
nents nonzero. The quasiquark spectrum is then easily 
determined analytically, and the dispersion relations of 
the individual fundamental colors, denoted as r,g,b, are 

Kr,g) ~ \/^k(r,g)> +^ ' k(b) = lsfc(&)l> 



E 



±, with 



£fe(r) =€h( g ) =£k + e[n 



Mg \ 



2^' 






■%£ _i_ /i-*-j^ 



(15) 



At zero temperature to which we will from now on limit 
our discussion, the density of an individual (fundamen- 
tal) quark color a is 



(«) 



r [ d A k e4fc (a 
la ~ 2-LJ (2tt) 3 El, , 

e=± J v ; k ( a ) 



(16) 



A. Calculation for fixed chemical potentials 



the momentum to zero, one finds 



With the particular orientation of the condensate that 
we chose, only the T§ generator of the color SU(3) devel- 
ops nonzero density, n% = (n r + n g — 2n(,)/(2yo). It can 
be made to vanish by tuning /itg appropriately. Let us 
assume that this has been done and calculate the propa- 
gator of the four modes that couple to Ai : 2- These form 
a complex doublet under the unbroken SU(2) subgroup 
and, had we not imposed the neutrality constraint, they 
would give rise to two type-II NG bosons with quadratic 
dispersion relation at low momentum 14, 16]. A neces- 
sary condition for these to appear is nonzero density of 
some of the color charges (see Schafer et al. in Ref. [La|). 
That is why one expects to recover four usual (type-I) 
NG bosons once the system is made neutral. 



Instead, we will show that once the color chemical po- 
tential is adjusted in order to make the system color neu- 
tral, the propagator of A 12 has a double pole at the fre- 
quency ojq — /i 8 y'3/2. This was in Ref. ]9] misinterpreted 
as a manifestation of two pseudo-NG states with degen- 
erate masses. On the contrary, it actually means that the 
system exhibits an instability of the same type as revealed 
in the scalar toy model in Section [II Al In order to under- 
stand this, one should note that the 2SC phase possesses 
an unbroken U(l)^ symmetry, corresponding to the blue 
quark number and generated by B = B — 2Tg/y/3 with B 
the baryon number. In the space of A a , which transforms 
in the antitriplet representation of SU(3), this symmetry 
is generated by the matrix diag(l,l,0). The two de- 
grees of freedom contained in the complex field Ai (and 
equivalently A 2 ) carry opposite charges, B = ±1; they 
are a particlc-antiparticle pair. Recalling the Umezawa- 
Kamefuchi-Kallen-Lehmann spectral representation, the 
particle pole should show up at positive frequencies, while 
the antiparticle one at negative frequencies. Since the 
chemical potential required to make the 2SC state neutral 
is typically negative, the double pole at w thus actually 
describes an antiparticle with mass |wo| & n d a particle 
with a negative mass, — |wo|- This is yet another man- 
ifestation of the seeming instability of the color-neutral 
2SC state with respect to fluctuations that generate off- 
diagonal color charges [6j. In the following, we will show 
how this problem can be fixed in a way that renders the 
theory stable and the NG bosons exactly massless. 



Let us now proceed to the proof of the existence of 
the double pole in the Ai propagator. This is defined as 
usual by U^(x - y) = -i(0|T{Ai(x)AJ(j/)}|0), where 
"T" denotes time ordering. The inverse propagator, or 
polarization function, IIn is most conveniently evaluated 
in the random phase approximation. Upon performing 
Fourier transformation to momentum space and setting 



iiiiH 



e=± 
E k{r) 



d 3 k 



1 



( 2 ^) 3 E Hr) 
e Sfe(r) 



(e£ W ) 



E 



fc(r) 



l£fc(6) 



(17) 



The polarization function II22 is identical. Using the gap 
equation at zero temperature, 






1 

e=± 



d 3 k 1 



(18) 



the inverse propagator (|17[) is easily brought to the form 

d 3 k 1 



IIii(w) 



0«(6)) 



(2^) 3 El (r) 



El 



IS 



El 



l& 



(19) 



J k(r) ' ISfe(6)l M fc(r) ls fc(fc) ' _ 

The prefactor in Eq. (|19p immediately tells us that 
there is a pole at w = ujq . Since this pole becomes exactly 
massless in the limit of fig — 0, thereby representing a 
NG boson, it is the particle pole. To find the antiparticle 
pole requires some further manipulation. Evaluating the 
expression in brackets at u = ujq yields 



s\/3 



El, N 

k(r) 

<?(«a ( 6)) 



< 



(b) 



2 



E] 



k ( r 



<l 



(i>) 



k(b)i 



(20) 



E k(r) 

E k(r) S S^« 



e *k(r) 

k(b)) 



E k{r) 



< 



k(r) 



< 



k(r) 



A 2 

Using Eq. (|16|) . one then arrives at the conclusion that 
the integral (including the sum over e) in Eq. (p~9|) equals 
— ns-\ /3/ (2A 2 ). In fact, it was already observed in 
Ref. [l7| that the coefficient of the term linear in to in 
the expansion of the inverse propagator around the par- 
ticle pole is proportional to the density n§. The present 
result is just a generalization to nonzero values of fi$ . We 
therefore conclude that when the chemical potential /zs is 
tuned so that n% = 0, the propagator indeed has a double 
pole at u — loq. 

We note in passing that, as the second line of Eq. (|2"0)) 
clearly shows, the double pole occurs outside the two- 
body continuum, so the implied instability cannot be al- 
leviated by decay processes. 



B. Induced fluctuations of chemical potentials 

We would now like to apply the same strategy as in 
Section in Al to remove the instability revealed above, and 



to show that all NG modes in the 2SC phase are exactly 
massless as they should. However, due to the complicated 
form of the mean-held thermodynamic potential in the 
NJL model, it is not possible to solve analytically for the 
color chemical potentials as a function of the (uniform) 
collective fields. Fortunately, this is not really necessary 
if we are only interested in the mass spectrum of the 
collective excitations. 

Let us consider generally a system possessing a set of 
(real) order parameters, A a , constrained by the require- 
ment of zero densities, n a , of a set of conserved charges. 
Introducing the associated chemical potentials, fi a , its 
thermodynamics is governed by the grand canonical po- 
tential, Q(A a ,n a ). The values of the order parameters 
and the chemical potentials in the ground state are deter- 
mined by the set of gap equations, dil/dA a = 0, together 
with the constraints, 



dHa 



= 0. 



(21) 



In order to study the fluctuations in a neutral system, 
one uses this equation to eliminate the chemical poten- 
tials in favor of the order parameters, and thereby arrives 
at a thermodynamic potential as a function of A a solely, 
£7(A a ,/i Q (A a )). Such a thermodynamic potential has a 
minimum that coincides with the simultaneous solution 
of the gap and constraint equations for Q(A a , /j, a ). In 
an unconstrained system, the mass matrix of the collec- 
tive modes is, up to an irrelevant factor, proportional to 
d 2 £l/dA a dA},. Once the constraint (|21[) is imposed, this 
must be obviously replaced with the total second deriva- 
tive, d 2 il/dA a dA D . Differentiating Eq. (f2~Tj) with respect 
to A a (this builds in the requirement that the constraint 
be satisfied for all values of the order parameter, at least 
in the vicinity of the equilibrium) , one arrives at 



d 2 n 



d 2 Q 



d 2 fl d\i a djip 



(22) 



dA a dA b 8A a dA b dfj, a dfxp dA a dA b 

d 2 n d 2 n (d 2 ny 1 d 2 n 

8A a dA b 8A a dfi a \dfidfij a/3 dA b dfifs' 

All partial derivatives are to be evaluated at the ground 
state values of the order parameters and chemical poten- 
tials. 

In the next subsection it will be demonstrated explic- 
itly that this prescription results in exactly massless NG 
bosons in the color-neutral 2SC phase. Here we just note 
that this is very natural: the potential Cl(A a , fi a ) is in- 
variant under simultaneous transformations of the order 
parameters and the chemical potentials. Therefore, the 
potential f2(A a ,/i a (A a )) is invariant under the symme- 
try transformation of A a , and the masslessness of NG 
bosons follows immediately from the Goldstone theorem. 

Before proceeding to the explicit NJL-model calcula- 
tion, we would like to point out one subtlety hidden in 
Eq. (22). The matrix d 2 n/d[i a dfi f3 in Eq. (22) repre- 
sents minus the density-density correlator (or the color 



number susceptibility matrix) in the ground state so that 
it is negative-semidcfinite. The charges of the unbroken 
symmetry have by definition sharp values in the ground 
state, and thus give rise to zero modes of the correlator. 
The matrix of second partial derivatives is therefore not 
invertible. There is a simple remedy: one adds a term 
— C 2 /i Q /i Q to the thermodynamic potential which makes 
all expressions well defined. In the end, the limit £ — > is 
performed. Since the unbroken generators do not couple 
to NG bosons, this subtlety does not affect the calcula- 
tion of their mass spectrum. 



C. Induced corrections to NG boson masses 

We will again focus on the complex doublet of NG 
bosons that couple to Ai, 2 [18j . Recall that these fields 
are complex and we are interested in the derivative 
d 2 Ct/dA\dAi which yields the static part of the inverse 
propagator (|T7|) . 

We need the three second partial derivatives of the 
thermodynamic potential, d 2 Sl/dA\dAi, d 2 ^l/d^ a d^p, 
and d 2 ft/dAld/j, a . The first one is already contained in 
Eq. (HI), 



1 o 2 n 

VdAldAi 



= -fj, S y/3X, 



X 



E 

e=± 



d 3 k 



sgn« 



fe(6)> 



(23) 
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The density-density correlator is given by two one-loop 
diagrams, with normal and anomalous components of the 
quark propagator. Carrying out the partial derivative of 
the thermodynamic potential (|13j) . we obtain 

d 2 U, 

Pj. q = Tr ("fQTaSqq'yoTpSqq — JoTaSqqJoTp Sqq) . 

The subscripts of the propagator denote matrix elements 
in the Nambu space. Inserting the quark propagator (T4) 
and carrying out the frequency integration, the first term 
evaluates to 



0V o 
^-A 2 tT(T a V 12 T V 12 )Y 

8V MT a V 12 T p V 3 ) + tr(T a V a TpVn)] Z, 
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where 
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and V12 and V3 are projectors on the subspace of first 
two colors and the third color, respectively. The sym- 
bol "tr" stands for a trace over color and flavor in- 
dices here. The second term in Eq. (|24l) reduces to 



-ifiV/2) tr(T a MTJ M^)Y , where M is a matrix in color- 
flavor space defined by M — 75-P3. Putting all the pieces 
together, we obtain 

1 d 2 n 



V d n a d up 

= -\ ( \ ^\ A 2 Y in the (T ,T 8 ) sector, (2?) 

= — Z5 a /3 for a = 4,5, 6, 7, 

and zero otherwise. Finally, the mixed partial derivative 
of $7 is given by a one-loop diagram with one normal and 
one anomalous propagator, 



d 2 n 



Tr (j T a S gg P a S g q) , (28) 



dA* a d^ a 
which, after appropriate manipulations, yields 

VdK^~ =tr(T Q M 75j P Q 7>3) X + ~ tr(T a M l5 P a V 12 ) 
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For the first component of the diquark field, the 
color-flavor traces are ti{T a M^^PiPi 2 ) = and 
tx{T a M^PxPz) = (0,0,0,-A,-iA,0,0,0) for a = 
1,...,8. 

Finally, we combine Eqs. {25} , {53}, {57}, and {55} to 
arrive at the simple result 

73^r-^i(»«->' x )- < 30 » 

A slight manipulation with the integrals reveals that the 
expression in the parentheses equals the charge density 
n%. This concludes the argument that once Eq. {22} is 
used instead of a simple second partial derivative to cal- 
culate the mass spectrum, the NG bosons are rendered 
exactly massless, as required by the Goldstone theorem. 



IV. DISPERSION RELATIONS OF NG BOSONS 
IN HIGH-DENSITY APPROXIMATION 

The approach taken in the previous section was 
straightforward, yet intrinsically limited to uniform or- 
der parameter fluctuations, so demonstrating the mass- 
lessness of the NG modes was the best it could do for us. 
Here we will follow the path outlined in Section IIIBI and 
construct the effective action for the NG modes in the 
2SC phase. We will use the high-density effective the- 
ory (HDET, see Ref. [1J| for a review), providing model- 
independent results in the limit of high baryon density. 

The HDET Lagrangian is constructed using quark 
fields and their coupling to colored diquarks. For the 



2SC superconductor, it takes the following form, 



A* a q^T 2 e a q + + A a q f + T 2 e a q*_ 
(L->R,A a ->-A„). 



(31) 



Here q stands for the left-handed positive-energy- 
projected quark field and the subscript ± denotes the 
(conjugated) velocity. The two four-velocities appear- 
ing in the Lagrangian are defined as V* — (l,v) and 
V^ = (1, —v); the sum is taken over pairs of patches on 
the Fermi surface characterized by velocities ±t>. Fur- 
thermore, t 2 is the Pauli matrix in the flavor space and 
\x is the chemical potential matrix in the color space; we 
consider the form p, = /il + Yla=i ^aT a ignoring the 
electric charge neutrality. Also e a — 2T^,2T^,2T 2 for 
a = 1, 2, 3 are the color antisymmetric generators, and ac- 
cordingly the A a 's stand for the background quark pair 
fields. The subscript a represents the color antitriplet 
indices which can be thought of as color-antisymmetric 
pairs, [g,b], [b,r], [r,g], respectively. Finally, the "L — > R" 
in the last line of Eq. (|31|) is to remind us that an equiv- 
alent Lagrangian for the right-handed quarks has to be 
added. 

By introducing the Nambu doublet field, ^> = 
(q + , —q*_) T , the Lagrangian {3T} becomes 



-HDET 



9^ \ -A*-nF„ iV-d-h* 



(L 



A* a r 2 e a iV-d-P 
R, A Q -> -Ac 



(32) 



Let us construct the effective Lagrangian for the NG 
modes. Following the logic explained in Section |H] we 
have to consider fluctuations of both the pairing field 
A a and the color chemical potential matrix jl, which is 
treated as a secondary order parameter induced by the 
color neutrality constraint. Without loss of generality, we 
adopt the convention that the ground state is character- 
ized by the nonzero condensate A3 = A3 = A. Accord- 
ingly, only the color chemical potential in the Ts-direction 
is nonzero in the ground state. The actual values of A 
and /is are set by the gap equation and the neutrality 
constraint. Having all this in mind, we parameterize the 
NG fields A a and chemical potentials /i Q in Eq. {35} as 



A a e a =U(7v)2AT 2 U(n) T , 
where U is expressed in terms of the NG fields tv as 



(33) 



U(n) 



exp 



1 



y^7T a (x)T Q 



Q=4 



(34) 



and we introduced the decay constant / w . In principle, 
an independent decay constant should be used for ev- 
ery real irreducible representation of the unbroken global 



symmetry, in this case one for a = 4, 5, 6, 7 and one for 
a = 8. We just use the same symbol for them in order 
to simplify the notation. 

It is now easy to obtain an effective action for the NG 
modes in a gradient expansion. Following essentially the 
same steps as in Section IlIIl we integrate out the quark 
fields to obtain 



S* cff =-Trlog5- 



(35) 



where S is the quark propagator in Nambu space in pres- 
ence of the fluctuating order parameters. This formula 
is to be contrasted to Eq. (flU)) . The mass term A a A*/G 
is missing here and consequently the NG nature of the 
order parameter fluctuations is made manifest. 

We find the following explicit expression for the prop- 
agator up to the second order in the NG fields, 



S~\x,y) = 



iV-d + fi -2Ar 2 T 2 
-2Ar 2 T 2 iV-d-p, 

J-rr 



5\x - y) 



(36) 



where 



-(i) - 
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— fisE a Ar 2 r a 

— M8S a/ 3 Ar 2 r a/ 3 



a/3 



(37) 



and S a , T Q , E Q/ 9, and z. a p 



are matrices in the funda- 
mental color space. The explicit forms of these matrices 
are found by expanding Eq. (|33|) up to second order in 
the NG fields. We note that the diagonal entries propor- 
tional to Us take into account the effect of fluctuations 
of the chemical potentials that was missed in previous 
analyses [9, Ha]. The gradient expansion of the action 
thus reads 



SW^-TrlogSo- 1 



x tr 



d4x n a (x)Mx) tf 
S {x,y)Y^S Q {y,x)T,f 



(2) 



S (x,x)E^ 



(38) 



with Sq L being the first term of Eq. (|36|). The first 
term above is proportional to the thermodynamic po- 
tential in the equilibrium. The second term gives the 
tadpole contribution stemming from the fact that we 
employed a non-linear parameterization of the collective 
fields. Omitting the part of zeroth order in the NG fields 
and working at nonzero temperature within the Matsub- 
ara formalism, the bilinear part of the effective action 



becomes, in momentum space, 



T 



Sett — — 



2 ^ 

N 



d 3 P 



ir a (-itt N ,-P) 



(2tt) 3 ' 
x IL a p(in N ,P)TT P {in N ,P), 



(39) 



where fijy = 2NttT and N is an integer. Reality of 
the NG fields in the coordinate space implies the con- 
straint on their Fourier components, n a (— ifijvi — P) = 
ir a (ittN, P). The inverse propagator of the NG fields is 
given by 



n aP (iQ Nl P)=—^ 



T ^ 
T 



tr 



d 3 q 

(2tt) 3 

d 3 q 



2So(iu n ,q)^ 



tr 



(2tt) 3 
T,^S (iu n +in N ,q 



So(iu n ,q) (40) 
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with w„ = (2n + 1)ttT being the fermionic Matsubara 
frequency. In the Nambu space, the fermion propagator 
takes the form 



of,-,, n\ - ( S q<i( iuj ni q) S qq {iu) n ,q) 
b {zu> n ,q)- l 5w ( iWn>g) S q q (i Un ,q) 



(41) 



where the individual elements are given explicitly by 



( r ,g) 



S q q(iuj n ,q) = (1 -B)S qq 
S qq (iu n ,q) = 2T 2 T 2 S<f g > ! >\ 
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q(b) 
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Lul + (v-q- fi r ) 2 + A 2 '' 
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(42) 



lUJ n 



(v ■ q - fi b ) z 

A 



q(r,g) 

qq ujl + (v-q- fi r ) 2 + A 2 ■ 

As before, B=\t- 2T 8 /y/3 = diag(0, 0, 1) is the pro- 
jector to the blue quark space, or in other words, the 
blue quark number in the fundamental color space. The 
chemical potentials for the red, green, and blue quarks 
are defined by fi r ^ g = fi + fis/(2\/3) and fib = H — fis/yS- 
Upon carrying out the trace in the color and flavor 
spaces and performing the Matsubara summation, H a /3 
turns out to have the following block-diagonal structure 
in the adjoint color space, dictated by the unbroken sym- 
metry, 



n, 



a/3 



TTCll 

11 tad 



n ch 
lAmix H 

n ch 



J-J-mix 

ch | TTch 
tad 



■IT 



(ff4,7r 5 ) 



TTCh 
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ch i TTch 
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n c 



(43) 



("•a, 17) 



(n^ ad + n n ) (ff8) . 



The superscripts "n" and "ch" distinguish modes that 
are respectively neutral and charged with respect to the 



unbroken quantum number B. In accord with what we 
already observed in Section IIII1 there should be four 
charged modes that form a complex doublet of the un- 
broken global color SU(2) symmetry, and one neutral 
NG mode. Using the unbroken symmetry, the propa- 
gator in the charged sector can be diagonalized in the 
basis (7T4 ± in?,)/ y2, {ttq ± iTTi)/y/2. The various con- 
tributions to the propagator are labelled by subscripts: 
"tad" refers to the tadpole contribution coming from the 
first term in Eq. (|4TJ|) , whereas "mix" refers to terms that 
mix different components ir a , thereby giving rise to mass 
splitting of modes with opposite values of B. As we will 
see below, the spacetime dependence of the propagator 
at long-wavelengths can be extracted analytically within 
HDET. 



A. NG mode in the neutral sector 

Let us first have a look at the neutral NG mode, 7r 8 , 
which corresponds to mere phase fluctuations of the order 
parameter and hence, based on the toy model analyzed 
in Section UH should not suffer from artifacts associated 
with the color neutrality constraint. Its dispersion re- 
lation can be extracted from the gradient expansion of 
n" ad + Xl n (iQ,N, -P). Explicit computation yields 



n- ad = — |A^(A, Mr ), 

6 J IT 



(44) 



with </> being the anomalous pair density defined by 
(<£?<#') = -(T 2 )ij(2T 2 ) ab (j>{A,fi r y, N { = 2 is the num- 
ber of flavors. In terms of the quark propagator, it is 
expressed as 



cf>(A^ r ) = -T^I-^S^(iu; n , q ) 



(45) 



On the other hand, the gradient expansion of II n (if2jv, P) 
becomes 



n n (iQ, N , p) = n n (o, o) - A(m 



N 



BP Z 



(46) 



with II n (0, 0) = — II" ad which guarantees that the mode 
is gapless in accord with the Goldstone theorem. The 
coefficients A and B are evaluated within HDET in the 
limit of vanishing temperature and at the leading order 

in A//i, ^ 8 //i as 



.4 



jVfjVo 



B 



N { N 



(47) 



where iVo = /j, 2 /2tt 2 is the density of states at the Fermi 
surface. In order for the NG mode effective action to take 
the canonical form with A — 1 , we set the decay constant 
to f% = NfN /3. The phase velocity of tt 8 becomes 



(48) 




which is the usual result reflecting the number of space 
dimensions [ffj . Since 7r 8 is a singlet under the unbroken 
global color SU(2) symmetry and carries zero B charge, 
it is not much affected by the color neutrality constraint. 
The above conclusion therefore remains correct irrespec- 
tive of the value of the background color charge density 
tt-8 as long as fis/fJ- is small. In the following, we will 
concentrate on the i?-charged sector of the NG effective 
action. 



B. NG modes in the charged sector 

We are now ready to analyze in detail the excitation 
spectrum of the charged NG modes in the (7r 4 , 7r 5 ) sector. 
[The spectrum in the (ttq, it-;) sector is identical as a con- 
sequence of the unbroken global color SU(2) symmetry] 
Before going into details, let us briefly pause to overview 
the general properties of the charged sector. Following 
the discussion below Eq. (|4"5)) . we define the combinations 
7r± = (-7T4 ± z7T5)/y / 2 which carry the charges B = ±1. 
The low-energy behavior of these can be extracted from 
the polarization function in the charged (n^, -^-sector 
in Eq. (|43|) . The components of the propagator sat- 
isfy the complex conjugation properties II ch (— iil^, P) = 
n ch (itt N ,P) and U m ^{-in N ,P) = -n mix (iQ N ,P). In 
the basis (ir + , 7r_ ) , the polarization matrix becomes di- 
agonal with the entries TI± = LT^ d + IT h qp «n m i x . The 
charge conjugation then implies that II_(iOjv, P) = 

n+Hfijv,P). 

In order to clarify the role of the fluctuations of the 
chemical potentials and the coupling between charge den- 
sity fluctuations and the NG modes, we decompose the 
functions IT h and n mix as 
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ch 
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IT 
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(49) 



where the first, second, and third terms are propor- 
tional to A 2 , A/is, and /x§, respectively. Diagrammat- 
ically, they are expressed as the quark one-loop diagrams 
containing two anomalous propagators, one normal and 
one anomalous propagator, and two normal propagators, 
is computed explicitly as 



TTCh 

il tad 



TTCh 

11 tad 



J n 



2 A0(A,/x r ) 



4/. 2 



2y u 8 n 8 (A,/i), 



(50) 



where n 8 (A, jl) = (n r + n ff — 2nh)/2y / 3 is the color charge 
density at an arbitrary /i 8 . In contrast to the tadpole con- 
tribution (|4"4")) in the neutral sector, there is an additional 
term proportional to /z 8 n 8 which comes from the tadpole 
diagram with a quark loop containing the normal propa- 
gator. This is because the charged NG mode is sensitive 
to the secondary order parameter, that is, /ig. This addi- 
tional term is absent if we ignore the fluctuations in color 
chemical potentials. 

We remark that 7r + ,7r_ correspond to A^A^ in the 
notation of Section IIIII Therefore, the curvature mass 



10 



matrix calculated in Section IIIICI can be obtained as the 
long-wavelength limit of the polarization function, 



1 d 2 n 

V dA* 1 dA 1 



lim n±(o,P) 

P— ¥0 



(51) 
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(o,o) + n ch (o,o 



We shall now calculate in detail the polarization func- 
tions II , and hence the dispersion relations of the tt± 
modes. In order to elucidate the role of the color neu- 
trality and the necessity to implement it carefully, we will 
analyze three different scenarios: 

(i) The case with fig — 0. In this case the primary 
order parameter, that is A, induce nonzero color 
charge density in the system. We expect the spec- 
trum to contain type-II NG bosons. 

(ii) The case of color neutrality ensured by "hard" 
(fixed) chemical potential(s). The chemical poten- 
tial fig is tuned so that there is no color charge in 
the ground state, but we ignore the chemical poten- 
tial fluctuations; this is analogous to the discussion 
in Section MI Al Technically this means that we 
discard the specific contributions to the polariza- 
tion function which couple to color density such as 
n^ 1 Ufl etc. In this case one may naively expect 
four NG bosons to acquire nonzero masses since the 
hard external background fig serves as an explicit 
symmetry breaking source in the quark sector [9[ . 

(hi) The full analysis of NG bosons in the neutral 
system. We take into account fluctuations in 
the chemical potentials and their coupling to NG 
bosons; this is analogous to the discussion in Sec- 
tion [nTCj Technically this means including all con- 
tributions to the polarization function. We ex- 
pect in this case to recover five massless type-I NG 
bosons. 



In order for the NG boson Lagrangian to have the canon- 
ical form, we adjust the decay constant as f 2 = NfN /2. 
The knowledge of the diagonal elements of the polar- 
ization matrix is not sufficient to determine the disper- 
sion relation or even to conclude that there are two mass- 
less modes. To that end, we need to evaluate the offdi- 
agonal component for which we obtain the result at the 
leading order in the gradient expansion, 



n m ; x (w,P) = -iCu, 



(54) 



where C = v / 3n 8 /(2f 2 ). In HDET, the color charge 
density takes the value n s = N { N a A 2 /(2 v / 3fi). This 
is suppressed by the small ratio A/fi and thus belongs 
to the next-to-leading order in HDET. Nevertheless, we 
keep the term C since the charge density rig is finite and 
can be evaluated without the high-density approxima- 
tion. The polarization matrix in the (774,775) sector then 
becomes, in the gradient expansion, 



-Alo 2 + BP 2 
+iCw 



—iCoj 
-Alo 2 + BP 2 



(55) 



Upon diagonalization in the (7r + ,7r_) basis, this leads to 
n±(w, P) = — Alo 2 + BP 2 =F Cw, resulting in the disper- 
sion relations 
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c J 



C 
A 
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(56) 



to second order in momentum. We can see that 7r+ , hav- 
ing the like charge as the medium, becomes a type-II 
NG boson with a quadratic dispersion relation, while 7r_ 
with the opposite charge acquires a gap proportional to 
the color density %. This is in agreement with the gen- 
eral theorems concerning NG bosons in systems lacking 
Lorentz invariance [1J, [l6l[l7|. Also, it is clear that the 
system is stable at least at the quadratic order in the 
derivative expansion. 



1. Analysis for the case (i) 

We set fig = and fi r = fi g = fi^ = \i. After ana- 
lytically continuing the polarization function to real fre- 
quencies, iClff — > w + id, and expanding in powers of u> 
and P, one finds, up to the second order, 



n ch (o;, P) = n ch (0, 0) - Alo 2 + BP 2 + • • • , (52) 

where II ch (0,O) = --i-A^(A,/i) which exactly cancels 

the tadpole contribution in Eq. (|50[) with fig — 0, re- 
flecting the Goldstone theorem. The coefficients A and 
B can be evaluated by employing the high-density ap- 
proximation (at zero temperature), which consists in the 

replacement J ,^ ? 3 — > N Q J^ d£ with £ = \q\ — /j, be- 
ing the momentum measured with respect to the Fermi 
surface. This yields the result |19j 
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(53) 



2. Analysis for the case (ii) 

The calculation is basically the same as in the case (i) , 
but we need to take into account the chemical potential 
differences, fi r = fi g 7^ fib- As before, we use the basis 
(tt + , 7T_), and thus need to analyze the polarization func- 
tion U+(lo,P) = Il£ d + W\lo,P) - iU mix (Lo,P). It is 
easy to evaluate the long- wavelength limit of II ch (lo, P) [= 

n c A h A (^,P)], 
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(57) 
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where £ = |g|— \i— /i 8 /(2y / 3) and £b — (■+\fi^a/2 are the 
momenta for red/green quarks and for blue quarks mea- 
sured from the Fermi surface. The first term in n AA (0,0) 
is cancelled by the tadpole contribution n^ d , but the sec- 
ond term survives. Therefore, in this case the curvature 
mass remains finite. 

Let us take a slightly different approach to the prob- 
lem. An explicit computation reveals that the func- 
tion H + (uj,P) depends on ui only through the combi- 
nation u + fib — Mr = w — ^3^/2- This fact sug- 
gests that it would be more natural to perform the ex- 
pansion about u = y3Ms/2 for H+(w,P) and similarly 
about uj = —y^3n$/2 for n_(w, P). In fact, we can eas- 
ily show that the offset of II ch (aj,P) =p iH m i x (uj,P) at 
uj = ±y3/Z8/2 completely cancels the tadpole contribu- 
tion so that 

11+ (w = VW2, o) = n (w = - VW2, o) = 0. (58) 

Therefore, we here perform the gradient expansion of 11+ 
about uj = \J?>pL%l c L. Up to second order in uj and P we 
find 



U+(uj,P)=A(uj- VW2) 

- B(uj - s/ln & /2f 



CP 2 



where A is proportional to the color density, 



A = 



A 



2 -n 8 (A,A) = 0. 



(59) 



(60) 



This is nothing but the neutrality condition which deter- 
mines the value of fig . In the high-density approximation, 
this can be calculated explicitly, /i 8 = —A 2 / (2y/3fi). 
(This is actually a next-to-leading order result, for it is 
suppressed as compared to A by the factor A//x.) One 
can also provide explicit integral formulas for B and C; 
since they are rather complicated, we again use the high- 
density approximation with the result 
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where N 0r = [i 2 ./2ir 2 . Setting f n = ^NqNi/2 and the 
phase velocity v^ = l/y3, the polarization functions 
acquire the final form 



IT 



-(w t73m 8 /2) 2 



(62) 



near uj = ±y3/xg/2, P = 0. 

We can see that under the color neutrality condition, 
the propagator of the charged NG modes acquires a dou- 
ble pole at uj = ±-\/3/j 8 /2 which clearly indicates some 
kind of instability as already explained in Section IHI Al 
These shifted poles were recognized as the "mass" of the 
pseudo-Goldstone modes due to the "explicit symmetry 



breaking" by /z 8 [HHl since detll(w, 0) ex (w 2 -3/4/4) 2 . 
Finally, we remark that, in accord with the discussion in 
Section Pi Al the instability is also clearly seen in the cur- 
vature mass squared, 



1 d 2 n 

uaA*aAT 



n±(0,0) = -%<0. (63) 



3. Analysis for the case (Hi) 

We finally consider the most general case with an ar- 
bitrary value of /k 8 and the corresponding color density 
n 8 , but with a proper account of the fluctuations in the 
color chemical potentials. 

Let us start with the diagonal element of the polariza- 
tion matrix, Lltad + II ch (u;,P). In this case we need to 
take into account the contributions from II A and EF~ 
in addition to n AA whose long-wavelength limit is given 
in Eq. (|57j) . Computing all the contributions using the 
explicit expressions for quark propagators, we obtain, in 
the long-wavelength limit 



n£ A (o,o) 



n 



ch 



„a(0,0) 

_7V f 
p 



IT 



2 A0(A, Mr ) 



c ^( > ) 
3 



VI 



2j u 8 n 8 (A,/i). 



(64) 



Note that the above offset is completely cancelled by the 
tadpole contribution in Eq. (|50|) . Thus, the diagonal en- 
try of the polarization matrix can be expanded up to 
second order in uj and P as 



iLad + n ch = -BUJ 2 



CP 2 



(65) 



where B and C are given by the same expression as in 
Eq. ([53)) up to the leading order in \i%j [i. This is true 
regardless of the value of /^ 8 as long as the color chemical 
potentials fluctuate with the primary order parameter 
according to Eq. (|3"3")l . In this case the curvature mass 
squared always vanishes, 



1 d 2 n 

V dAJUAT 



= n tad (o,o) + n ch (o,o) = o 



(66) 



In order to understand the low-energy behavior of the 
NG modes, we still need to evaluate the off-diagonal part 
of the polarization matrix, n m j X ; this consists of three 



parts LT 



AA 
mix' 



n 



Aii 



and n^ x . 



Each of the contributions 



takes quite a complicated form, but surprisingly putting 
them all together gives rise to the following simple for- 
mula at the lowest nontrivial order in the gradient ex- 
pansion, 



T M A 



n££( W| P) + C?>, P) + rO, P) 

,V3 



~i^-ujn s (A,fi)+0(uj i ,ujP z ). (67) 
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Lagrangian is to introduce two matrix-valued fields 



(n 8 > 0) 



(n 8 < 0) 



FIG. 1. (color online) . The schematic picture of the NG boson 
masses as a function of the color chemical potential /is . The 
neutral NG mode is always type-I and it is always gapless to 
the lowest order of approximation. 



This only depends on fig through the color density ng. To 
summarize, the polarization matrix in the charged sector 
behaves at long wavelength as 




CP 2 



-iYfrn s (A,fj,)u}' 
-Buj 2 + CP 2 



(68) 



It is now obvious that only when fig is tuned such that 
the color density ng vanishes, the full set of five type-I 
NG bosons are recovered. 

In Fig. Q] we show the schematic plot of the excita- 
tion gaps (masses) of the NG modes as a function of fig. 
There is one neutral NG boson which is always type-I ir- 
respective of the value of fig and the charge density of the 
system. The four charged modes with B = ±1 are, on the 
contrary, sensitive to fig. When fig > — A 2 /(2-\/3/z), the 
system is positively charged with color ng and the two 
type-II NG modes appear in the B = +1 sector; their an- 
tiparticles have the gap u> — y/3ng/(NfN ). On the other 
hand, if fig < —A 2 / (2-^/3^), the system is negatively 
charged, and the quantum numbers of type-II NG modes 
and those of their massive partners are interchanged. At 
the neutrality point, fig — —A 2 /(2^fi), the type-II NG 
modes change smoothly to the type-I NG modes. 



C. Effective Lagrangian for the charged NG modes 

The results for the dispersion relations of the charged 
NG modes can be conveniently encoded in an effec- 
tive Lagrangian. Generally, the nonlinear effective La- 
grangian can be constructed using the coset field U(tt) 
introduced in Eq. (34). A symmetry transformation 
g S SU(3) x U(l) acts on this field as U(ir) A U(ir') = 
gtW(7r)/i _1 (7r, g), where h(n,g) is a suitable element of 
the unbroken subgroup SU(2) x U(l)^. 

A convenient way to ensure invariance of the effective 



T = UT 2 W 



A = UTdA\ 



(69) 



These embody the existence of two order parameters, 
transforming in the antisymmetric rank-2 tensor repre- 
sentation and the adjoint representation, respectively. In 
the context of HDET, they were introduced in Eq. (|33l) . 
see also Eq. (TT0|) for the analogous construction within 
the toy model of Section [TTJ Under the action of g, 
these fields transform homogeneously as T — > gl~g T and 

A A gAg*. 

As long as we are interested only in the dispersion re- 
lations of the charged NG modes, we can use any of the 
fields T, A to construct the effective Lagrangian. (They 
lead to kinetic terms which are equivalent up to a re- 
definition of the parameters of the Lagrangian.) We 
choose A without loss of generality, and write down the 
most general SU(3) x U(l)-invariant Euclidean effective 
Lagrangian up to second order in the fields, 



C, 



ff',0 



4/2 



tr (d T Ad T A + vlVA ■ VA) , 



(70) 



where the prefactor has been chosen just for convenience. 
Since there is an unbroken U(l)g symmetry, we should 
allow for the possibility that it is endowed with a chemical 
potential. This is done by the replacement 



d T A^d T A- fi s [B,A}. 



(71) 



Finally, the background ng, charge density breaks the 
global color symmetry "explicitly" . Since the charge den- 
sity transforms in the adjoint representation, we have to 
add a term 



SC C 



--M 2 f 2 tr(AT s 



(72) 



where the parameter M 2 measures the amount of explicit 
symmetry breaking, very much like the pion mass in the 
chiral perturbation theory of QCD. Expanding the action 
up to second order in tt and discarding irrelevant terms, 
we arrive at the lowest-order action for the NG modes in 
the (774,775) sector, 



C 



off 



= tr [(9 t tt - ^^ [B, 7v}) 2 + vI(Vtv) 2 + M 2 tt 2 
--[{d T + n s )ir-][(d T - Mb)tt+] + ^ Vtt_ • Vvr 4 

+ Af 2 7T_7T + , 



(73) 



where now 7r = 7T4T4 + ir^Tc,. The low-energy couplings 
fig, v 2 , and M 2 have to be computed in some micro- 
scopic model, as we have done above using HDET. Let 
us emphasize that fig is to be interpreted as an effective 
chemical potential; it is not equal to the chemical po- 
tential of the blue quarks even though B represents the 
operator of blue quark number up to an overall factor. 

The Lagrangian (|73p implies the dispersion relations 
for the 7r± modes, 



^/vIP 2 +M 2 t^b- 



(74) 



Concretely, in the case (ii) where the fluctuations of 
color chemical potentials are ignored, we obtained /i^ = 
-y/Sfjis/2, vl = 1/3, and M| = 0, see Eq. (JB2J). The in- 
stability of the 7r + mode, akin to the instability revealed 
in Section [II A[ is now manifest. 

The cases (i) and (hi) can be treated together; the 
expressions valid for the former can be obtained by set- 
ting ^ 8 = in those for the latter. We obtained fi B = 
V / 3n 8 (A, / i)/(4/2) = y3 n8 ( A ,A)/(2AWo), t£ = 1/3, 
and M„ = \n B \, see Eqs. (|55|) and (|58|) . In the case (i) 
where color neutrality is not imposed, both \ig and M^ 
are nonzero and we find one type-II NG mode and one 
massive mode with the gap 2M n . On the other hand, in 
the case (iii) both parameters are zero and we find two 
type-I NG bosons with the phase velocity v n . 

On the general note, the above result can be inter- 
preted as the background charge density in the fermion 
sector acting as a chemical potential for the modified 
baryon number. This provides some intuitive under- 
standing of how the propagation of the NG modes is af- 
fected by the background charge. Putting together the 
collective modes in all (charged as well as neutral) sec- 
tors, the conclusion for the two physical cases (i) and 
(iii) is as follows. For n§ ^ 0, there are three NG modes, 
two of which are type-II with B = 1 and one type-I with 
B = 0. When n§ = 0, two gapped modes with B = — 1, 
the antiparticles of the B — 1 NG bosons, become gap- 
less, and five type-I NG modes are correctly recovered. 
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well as explicit breaking of the global color symmetry re- 
ported previously in literature disappear. To avoid confu- 
sion, we first addressed the problem in Section [TTT] using 
the same setting as in the previous publications. Only 
then, in Section IIV1 we provided a mo del- independent 
approach to the dispersion relations of the NG modes. 
While we used the 2SC phase as a specific example, it 
is obvious that our conclusions are valid more generally. 
Apart from other color-superconducting phases, they ap- 
ply equally well to all systems with a set of conserved 
charges, some of them being demanded to be zero in the 
equilibrium. 

It should be pointed out that the problem investigated 
in this paper is to some extent artificial. In the full QCD, 
color neutrality is automatically guaranteed by the gauge 
dynamics. Moreover, the collective modes discussed here 
are all colored and thus do not correspond to physical 
states in the spectrum; they will eventually be absorbed 
by the gluons by means of the Higgs-Anderson mecha- 
nism. It was shown that the collective excitations, the 
NG modes, play an important role in establishing gauge 
invariance of the Meissner effect [20] - We thus believe 
that our work provides a first step towards the proof of 
gauge invariance of the color Meissner effect in color su- 
perconductors, which has not been demonstrated satis- 
factorily so far. We leave a deeper consideration of this 
issue to future work. 
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